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CONTROL OF SOLID-APOGEE MOTOR DUPING ORBIT TRANSFER MANEUVER 


Zhang Yuntong 

(Beijing Institute of Spacecraft System Engineering) 


Abstract 

This paper provides a method for determining the optimum time and 
attitude in process of geostationary satellite launch at which the 
solid-apoaee motor to be fired. 

Subject Terms: firina time, satellite attitude, optimum design. 

Apoaee enoine 

(The manuscriDt was received on July 13. 1989.) 

I. Forward 

In the process of launching Geostationary satellites, the Apogee 
enoine on the satellite provides an increment in the velocity of the 
satellite which, ultimately, transforms the orbit the satellite 
travels from a transitional orbit into a geostationary orbit. The 
direction of spin of a geostationary satellite is the direction of the 
velocity increment. Therefore, the orbit transition control of an 
Apogee engine is to determine the firing time and the Spin direction 
(the so called firing attitude) for a given transitional orbit. 

Different firina times and firing attitudes can result in different 
orbits after 'rhe engine is stopped. Hence, the firing time and firing 
attitude should be carefully chosen to obtain the optimum orbit. In 
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the actual engineering practices, the choice is always subjected to 
some limitations, and. therefore, this problem can be thought of as an 
constrained optimization problem. In this paper, we will introduce: 

(1) the working model of Apogee engines, (2) the definition of 
optimization target functions, and (3) method of optimization. 

II. The Working Model of Apogee Engines 

1. Equation of Motion 

The equation of motion of the center-of-mass of a satellite during 
the working period of Apogee engine can be simplified as: 
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where x. y. z are the coordinates of the satellite in the earth 
equatorial coordinate system; ¥ is the distance from the satellite to 
earth's center: * is earth's gravitational constant; F is the thrust 

of the Apogee engine; m^ and m are the initial mass of the •stellite 

’ • -£• : ; 

and the rate of mass flow of propellent (per second); fg|i, «e ; Imd are 

the longitude and altitude of the firing attitude. The 4aMMNp*St£'on 

* 

should be carried out in the time interval of the Working period of 
the Apogee engine; namely, the integration starts at the teajg-.efhen the 
Apogee engine is ignited and stops when the Apogee engine is stopped. 
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After the transformed orbit is selected for a firing time t the 

initial conditions for the set of equations (1) can be obtained and F. 
mQ. m and the working time of the Apogee engine are all fixed 
oarameters. If a set of.* and i were chosen, the orbit of the 
satellite after the engine is stopped can be obtained from the 
integration of equation (1). In other words, for a set of firing time 
and firing attitude (tg the orbit of the satellite after the 

engine is stopped can be calculated. 

2. Pulsed Orbit-Transition Model 

For solid-Apogee motors, the working time is on the order of a few 
tens of seconds. During *its working period, the satellite is in the 
vicinity of the Apogee point. It can be assumed that the satellite 
obtains the velocity increment instantly. The position and velocity of 
the satellite at ignition and stop of the Apogee motor should be 
expressed in terms of the subscripts "0" and "f ". Their relationship 
can be simplified as 

*. = */ 

».=*/ 

*.=»«/ 

g. * V omm a a at 

f,=f,+Voemdmna 

i. = z % +vm8 

where V is the velocity increment provided by the Apogee motor and it 
guantity is: i-.c 

_ *. k . * 

where gg is the gravitational acceleration on the surface of the 
e»rth: I S p is the specific impact in vacuum; and T is the working time 
of the motor. Similarly, if the firing time tg is determined. x 0 . yg, 
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zq. xq. yo > ZQ can be calculated from orbit measurement and a and 6 
are the quantities which should be chosen. The relationship between 
the parameters x s . y s , z s . x s , y s . z s and the orbital parameters is 
determined bv the well known equations as those in [1], 

3. Comparison and Conclusion of The Two Models 

In the two methods mentioned above, method (1) is closer to the 
actual motion. When the working time of solid motors is a few tens of 
seconds, the difference of the calculated target functions can be 
ignored. The minimum difference in the target function is obtained 
when the difference in the firing times based on these two models is 
0.55 times of the working time of the motor. Therefore, optimal 
condition can be obtained if the actual firing time is 0.55T in 
advance of the time based on the pulsed model. 

III. Taroet Function 

1. Definition of the Target Function 

In order to evaluate the advantages and disadvantages of the chosen 
parameters tg. a , and 6, the optimal target function is defined as: 
the corrective velocity increment necessary to bring the orbit at 
engine stop to the geostationary ncbit. The optimum set;of functional 
variables is the set of tg. a, and <5 uhieh minimizes the target 
function. 

2. Calculation of the Target Function 
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The orbital parameters at engine-stop are: geometrical longitude a* 
engine-stop A,, halt lono-axis of the orbit *.» orbital eccentric 
ratio e^: oibital tilt angle equatorial longitude of the up-ncmo 

point Q,. The parameters of the geostationary orbit are: fixed 

position geometrical longitude half long-axis of the orbit «* : 

orbital tilt angle i Q : equatorial longitude of the up-rising poim 

In order to maintain the north-south position should be around 

» - 

270° and i c should be about 0.1°. To maintain the east-west position. 
e s should be less than a certain value and this value was taken as 
zero when the target function was calculated. The corrective velocity 
in the target, function was expressed in terms of three parts: 

(1) Corrective velocity AK, necessary for the correction of half 
lone-axis and eccentric ratio: 

To correct e s to •..(eg-0), a dual-pulse was required and the 

veiocitv was increased (or decreased) at the hunch point. The 
ca 1 cu 1 at ion* 1 eauatmn? are: 


r. r = «.(!-«,) 



The corrective ve 1 ocity VJi ♦ F» — VJ 

(2) Corrective velocity AV2 Tiecessary for fixed-point capture: 
Assuming that 4 »> 4 * and in order to maintain the satellite at 
position a % . the satellite should drift towards the west after the 
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motor wjf stopped. If the rated drifting velocity is i (degrees/day), 
1<0 . The half long-axis corresponding to the drift should be 





The calculational equations for corrective velocity 
various and e s are tabulated in table 1 : 


Table 1: Calculational equation for corrective 

•1 1 




4*. 


«.>a< 


a,Cat r„>a„ ] 


r .p<-0o 1 «[>0j 


»>z (a<-a.) 
0 






ikev 1 : condition 


y 

It: table 1. an = f is the anaular soeed of earth's 

t>c 

(3) Corrective velocity At, necessary for orbital 
Orbital plane capture is achieved by correcting ig 
and the orbital plane should be tilted angle. 

The corrective veloc ity can be appro****** •«: 
when ip < t/ 2 . AV 3 -V()ip 

when i p > w / 2 . ^ 3 -Vo (w-ip ) (i p in radians) 


kv t according to 

velocities 


spin and 

plane capture: 

•A 4 to i 0 and Q% 1 

■' * 


If A s tA q . similar results can also be obtained. 





3. Target Functional Values 

The target function is A V- AVj+ Al^^-AVg and the relationship 
between the chosen parameters (t(j. a ) and the target function is 
explained in the above sections. 

IV. Optimization Process 

1. Method of Optimization 

Under the normal conditions, in the process of satellite iaunchina 
the error encountered m the transition orbit is usualiv very small 
and the reaion for optimization is also a very narrow reaion: on th*- 
other hand, the calculation based on the pulsed model is also very 
simple. Based on these reasons, the method of optimization can be the 
reliable listincr method and the searchina can be carried out in the 
vicinity of the standard values. 

2. Calculation of the Standard Values 

It was assumed that the standard values satisfy the two conditions 

<♦ 

1) i s> fi g equal ig, fig. respectively; 2) And V g are mutually 
perpendicular and x s x s +y s y s +z g z s -0. Further, it was assumed that the 
instantaneous roots at the uD-rising point before ignitloh^f the 
motor were A3. 63, 13, 103, 0 3, K3, fSl^lacing the instantaneous roots 
of at the ignition point. The calculation of the standard values of 
the firing time and firing attitude was then based on the familiar 
two-body equations: 
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006), = oo*r 0 oo8«, + sun,«ni.oo»<C,-a,) 

«n#, e = an (£? s — Q») tint,/tint, 

(y * I "**)-© 

Af = Afj +Et — f s sin Ei 


t =M 
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«’»*=P»/<l + r»C 0 B/,) 



e,= «m-‘Cv/ MP~»/r»V^ (f*>* # £» ** 11 * 111 ) *<S> 
<5» = sin -l ( sin*,sinC«i» + £j)) 
«,=C,+tg' l (oosi 1 tg(ji M +e I )) (tg -1 $)*RI-5 
008 9 % — sine, cos*/ 
oos r — - ctgfjctgy. 


sin ^4 = ^ «n P, 

' 4 

C0S>7= —tg3jctg(«j 0 + E») 

an 6 = sm5jCos(9 r , + 9 r ,) + oosA,sin(l r , + F4)co8(?7+x) 
sin(dj- a ) = sin (*%♦*«) sin (if+*)/oos A 


key - 1: £ 3/2 and E 3/2 are in the same quadrant. 

2: When f 3 is greater than »;.e 3 lies in the second quadrant. 
3: tg -1 and U 3 Q + e 3 lie in the NM quadrant. 

3. Conclusion of the Optimization Process 

The control of the orbit transfer maneuver of Apogee motor can be 
concluded as: 

( 1 ) measurement of the transition orbit: 
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(2) calculation of the standard values '’7, fc, 7, 

<-confirmation of the region-of-choice and step-length of t.o . 
and 6 : 

(4) calculation of the parameters a— . e 8> i # . after stop of the 
motor for each set of t. a^and i ; 

(5) calculation of the target function A V; 

(6) repeat of steps (3) «-*(5) to determine the t, a , and a 
corresponding to the minimum A V; 

(7) calculation of t-0.55T and (a.. & ) for actual firing time and 
firing attitude, respectively. 

* 
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APPLICATION 


OF SPLIT-PHA C; E CODE P9K SIGNAL IN DATA TRANSMISSION OF 

METEOROLOGICAL SATELLITE REMOTE SENSING IMAGE 

Dai Ying 

(Shanghai Institute of Spacecraft Engineering) 


Abst ract 

In this paper, split-phase code PSK signal is analyzed, its 
character is reviewed, the basic formation method is briefly 
described. The application of the signal in digital image transmission 
of astronautica1 remote sensing is introduced. And finally, photos of 
real-time data transmission are presented. 

Subject Terms: phase shift keying, phase modulator, picture 

transmission, meteorological satellite. aDplication 

(The manuscriDt was received on May 5. 1989.) 

(The definition of split-phase code can be found in reference 1.) 

I. Split-Phase Code Phase-Shift Keying Signal 

Split-phase code phase-shift keying (PSK) signal is the modulated 
signal of the carrier wave, based on the split-phase code signal, 
after phase-shift keying. It is a binary modulated signal; namely, the 
PCM-PSK/digital split-phase, which is also a primary phase-modulated 
signal without the auxiliary carrier waves. Its mathematical formula 

1 S : 

»()«(.« (1) 
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where u m - amplitude.- 

&>,- - frequency of the carrier wave ; 

Ap- - split-phase code phase-shift amplitude. 

c(t) - split-phase code, its value should be taken as +1 or -1 
randomly. 

The relationship between the amplitude and phase-shift angle of the 
split-phase code PSK signal is: 

yiS' ( 2 ) 

where u cm is the voltage amplitude of the modulated carrier wave. 

In the process of product research and development, according to 
equation (2). the phase-shift amplitude of the carrier wave by the 
random solit-phase code can be measured with the frequency 
spectroscope. Using this method, the modulation of the transmitter can 
be measured with hiaher accuracy. This equation can also be used in 
the orbital test of satellites. 

The mathematical expression of the frequency spectrum of split- 
phase code PSK signal is: 

F(a> -*) c ) = jo * 1 A? y *•) 



- r , r .... ■ 


where $ is a function with screening characteristics which is also 
called the 6 (function and T is the period of the split-phase code. 

According to equation (3), the frequency spectrum can be plotted as 
shown in figure 1. The shape of coverage of the spectrum of this code 
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is similar to that of the fundamental band signal frequency spectrum 
of the split-phase code. Therefore, in order to avoid interference, it 
was only necessary to filter the visual wave, before modulation, to 
accomplish the intra-band constraint of the spectrum of split-phase 
code PSK signal. Technically, this method is easier to perform than 
the intra-band constraint of the microwave after modulation and. also, 
this method is more economically feasible. 

The split-phase code PSK signal has certain characteristics which 
makes it very powerful in the digital image transmission. Firstly, 
since the phase shift anale of the split-phase code was not equal to 
it/ 2, the carrier wave part could be saved after modulation. Secondly, 
since the part of the fundamental frequency spectrum of the split- 
phase code at zero freauenev was zero and the phase-shift modulation 
of the carrier wave was the linear transformation, the part of the 
carrier wave after modulation existed independently and the 
interference of the chancre of the messaoe spectrum did not exist. This 
led to the features of stable carrier frequency part and the hiah 
signal to noise ratio. For example: when the narrow band in the 
vicinity of the carrier frequency in the frequency spectrum of the 
split-phase code A/=o.S5jHtf, the power of the message part of the 
continuous spectrum was less than 0.03* of the total power of the 
message spectrum and the actual signal to noise ratio could be as high 
as 75dB. 

A ' * 

In figure 1. the percentages were the results of the numerical 
integration of equation (3) and indicated the power of each specific 
wave peak as compared to the total power of all the wave peaks in the 
continuous wave spectrum. One can see that most of the power of the 
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1: Frequency spectrum of the split-phase code F'SK siqnai 
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continuous wave spectrum existed in the first wave peak which amounted 
to 90%. 

On the other hand, under the condition of total symmetry of the 
split-phase code, there was no dispersion frequencies of the modulated 
frequency spectrum of the split-phase code PSK signal. Also, there was 
no even parts. In this way. the erroneous tracking of the carrier wave 
by the ground station could be avoided. 

With these special features, the split-phase code P9K signal was a 
useful tool in determining the modulation functions of a transmitter 
and the design of the match of system frequency band. If the band 
width which occupied 90% of the energy was 2.8MHz, considering 1) the 
frequency deviations caused by the carrier wave and the signal itself, 
2) the acceptable signal to noise ratio and, 3) the deformation of the 
waveform of the signal, the band width of the receiver was chosen to 
be 3MHz . Even though the edge of this band width overlapped with the 
corresponding frequencv spectrum of the TIROS-N satellite of USA, no 
interference should occur because of the separation in space and time. 

If the special features of signal mentioned above were linked with 
the system function, the advantage of the split-phase code PSK signal 
was made even more obvious. For example, to transmit digital image 
with split-phase code PSK signal, when the receiving/demodulat ion 
equipment received the PSK signal, the PSK signal power which entered 
into the carrier wave lock-up phase loop of the receiver was very 

lb;- * • 

limited. In other words, even though the signal received flips the 
modulated signal, the signal to noise ratio of the part of the carrier 
wave which entered into the loop was high. In this way, the signal 
could be transmitted under the condition without the modulation of 






auxiliarv carrier wave and the receiver could perform tracking of the 
carrier wave and de-modulation of the sianal simultaneously. 

The disadvantages of the split-phase code PSK signal was that when 
the signal was transmitted through wireless transmission, the 
interference would occur due to the effect of multi-path electrical 
transmission. Therefore, when the system was tested or measured, the 
geographic effect should be adequately taken into account. 

II. The Formation of the Split-Phase Code PSK Signal 

The split-phase code PSK signal was iormed oy Keying the phase of 
the carrier wave signal with the split-phase code. In the following, 
two basic methods for the formation of split-phase code PSK signal 
would be bneflv explained and mathematical methods were used to 
analyze the formation process. 

Method 1: 

The basic principle is shown in figure 2 and the relationship 
between various phases is shown in figure 3. 

Fiaure 2 shows that the multiplier does not shift the phase of the 
carrier wave signal u(tl in the first half of the cycle of the split- 
phase code operation time, but in the second half of the operating 
cycle, u (t) is shifted 180 degrees. The mathematical analysis is: Let 

r 

u ( t ) “Umcosttot, when the output signal of the-aultipiier is 

» 

where k is a constant. 

At the same time, the output signal of the adder is 
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Fiaure 2: solit-phase code phase shifting principle 
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Fiaure 3: Phase relationship of various components 
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(A) 


S(») *».«*•*» 

A <p) 


where 


A m = y/ni+*!. t 
A* =aretg |^|. 


Therefore, the phase-shifting modulation of the carrier wave by the 
positive electrical potential (subsidiary code "l") of the split-phase 
code vos accomplished. When y<»<r/ u(t) was shifted 180 degrees. 

«l (#)**».«•<•*#+ -y+jr) 

s(t) = u«) + a,(Q 

( 5 ) 

= ».00M%f+it,*4natt 
= i4.0M(a>^- A?) 


where 


A.= >/«! + «}. > 


Af>= arctg 


tii 

a. 


Therefore, the phase-shifting modulation of the carrier wave by the 
negative electrical potential (subsidiary code "0") of the split-phase 
code was accomplished. It can also be seen that by controlling the 
accuracy in phase-shifting and the amplitude ratio of the two signals, 
the split-phase code PSK signal can be formed. 


Method 2: 

Figure 4 shows a reflective-type modulator which accomplishes the 
phase keying control of the carrier wave by controlling the path L of 
the carrier wave. When the split-phase code was at the positive 
electrical potential, the PIN tube was connected and point A was 
arounded and the carrier path was Lj . When the split-phase code was at 
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Fiaure 4: A reflective phase modulator 
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the negative electrical potential, the PIN tube was disconnected and 

the path of the carrier wave was L„. If Afi (L is satisfied, 

4 lei 

the carrier phase transformed fromf>,to fit j can be obtained. was the 
initial value. 

In actual application, the 3dB bridge was used in replace of the 
cyclotron in figure 4. In fact, the modulator manufactured in this way 
would perform better and actual manufacturing was accomplished. 

III. The Application of the Split-Phase Code PSK Signal 

The split-phase code PSK signal was used not only in the digital 
communication of spacecraft, it was also used in the transmission of 
high resolution remote-sensing images of TIROS-N meteoro 1 og 1 ca1 
satellite. The met eoro 1oalca1 satellites manufactured in China also 
used the split-phase code PSK sianal to transmit ultra-high resol’d ion 
cloud maps. This system consisted of two parts (satellite transmission 
and oround ie reiving > : 

Satellite Transmission: The satellite transmission included the 
electronic circuitries of the remote sensing system and the 
transmitter and antenna of the transmitting system. The remote-sensing 
system installed on the satellite was a five-channel scanning 
radiation monitor which received the thermal radiation data from the 
emission and reflection of the earth-atmosphere system and the visible 
light data which included the meteorological data from cloud top, 
earth surface and ocean. 

The remote-sensing equipment first transformed the gathered data 
into electronic signals and then quantized and digitally processed 
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these signals. After coupling the service code, the information was 
edited and according to the special format requirement, such as the 
service code of synchronous signal, satellite attitude, earth 
meteorological data, and the error message detecting sequence, the 
binary code was transformed into the split-phase code. After phase- 
shifting keying of the carrier frequency of the transmitted cloud map 
data, the split-phase code PSK signal was generated. This signal was 
transmitted co ground, via antenna, after amplification. 

Ground Receiving: After the ground station received the split-phase 
code PSK signal, the high-frequency part was amplified, the 
frequencies we>e mixed, intermediately amplified, and then transmitted 
to the lock-phas® de-modulator. During de-modu1 ation. the phase of the 
siunal was compared with the phase of the split-phase code PSK sianal 
carrier wave reference signal and the split-phase code signal was 
obtained. The amplified split-phase code was transmitted to the bit 
svnchrom per to extract the synchronous clock signal and the split- 
phase '-ie was transformed into the non-zero code and transmitted to 
the frame synchronizer. In the frame synchronizer, the beginning 
of a frame sianal--the 60 bit pseudo-random signal was examined to 
extract the frame synchronous timing signal and the cloud maps or the 
earth meteorological signal was separated into five channels. For the 
error message signal at the end of a frame of signals, the error rate 
would be examined and rate of error message of the whole string of 
data would be derived. If the error rate satisfied the requirement, 
the cloud map signal would be transmitted to the 4ata processing 
system for real-time image processing. Or, the data could be stored 
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for future processing to obtain the meteorological data of the cioug 
top. e^rth surface and ocean. 

Figure? b and 6 show two processed high resolution cloud maps which 

were transmitted by the high resolution cloud map transmission system 

(HRPT system' of the "Fong-Yun" No. 1 meteorological satellite via 

split-phase code PSK signal. Figure 5 shows the cloud map transmitted 

bv the HRPT system at 08GMT on September 8. 1988. It shows the rain- 

producing cloud system of the south eastern part of China. In the 

southern part of China, there were several cloud strings which fed to 

the rain-producing cloud region. The dark region at the bottom of the 

fInure was the South Sea and Hainan Island was at the center of it. 

Figure 6 shew? the visible light high resolution meteorological cloud 

mer* transmitted bv HRPT svstem. at the downstream of the Yangtze 

river at 6:30 GMT on September 20. 1988. This map shows clearly the 

-! ovd •list rib'it ion and the sand content of the downstream or th* 

Yangtze rive and the East Sea. It also shows the position of the five 

uiu* ; . *• tp gi’ the w i d t h of each string at the coast a 1 region vf t he 

CianoSu province. Also, the south pointed dispersion of the fresh 

water from the Yangtze river in the ocean can be clearlv seen. This 

dispersion of the fresh water affected the HangJu bay. Further, the 

mixing and distribution of the fresh water and sea water can be seen. 

(Figures 5 and 6 are in the third cover page.) 

All of these figures have clear images and clear geographic 

boundaries. The resolution was less than 1km. One of the high 

r 

resolution cloud maos was used by the "AVIApN WEEK" magazine. 

Data transmission was one of the key issues to produce high 

auality cloud maps. Therefore, it should be noted that transmission of 
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the remote-sensincr images by the split-phase code PSK signal was 
successful. These cloud maps provide important data for weather 
forecasting, typhoon monitoring, ocean detecting, and agricultural 
evaluation. 
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